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SECTION  1 


INTRODUCTION 


1-1  BACKGROUND. 

High  energy  neutrons  can  cause  upsets  in  modern  complex  microcir¬ 
cuits  at  exposure  levels  much  less  than  necessary  for  significant  permanent 
damage.  Electronic  systems  using  microcircuits  are  hardened  to  the  tran¬ 
sient  effects  associated  with  the  pulsed  ionization  and  permanent  damage 
effects  of  the  total  ionization  and  neutron  fluence  of  the  nuclear  weapon 
radiation  environment.  The  question  is  whether  current  hardening  tech¬ 
niques  are  sufficient  to  protect  a  system  from  deleterious  single  event 
upsets. 

Recent  experiments*^  in  Single  Event  Upset  (SEU)  phenomena  have  ver¬ 
ified  the  susceptibility  of  nMOS  dynamic  RAMs  and  nMOS  microprocessors  for 
high  energy  (nominally  14  MeV)  neutrons.  The  mechanism  is  ionization  in  a 
critical  volume  of  a  memory  cell  deposited  by  an  alpha  particle  from  an 
inelastic  collision  between  an  incident  neutron  and  a  silicon  atom.  The 
(n,a)  cross  section  is  -100  mb  for  neutron  energies  above  5-6  MeV.  A  14 
MeV  neutron  can  produce  alphas  with  energies  between  7.8  and  11  MeV.  Only 
2.27  MeV  of  deposited  ionization  energy  in  silicon  is  required  per  0.1  pC 
of  charge.  The  critical  charge  for  SEUs  in  many  VLSI  devices,  both 
memories  and  random  logic,  is  in  the  range  of  0.1  to  1  pC6*?. 

Since  it  has  become  well  established  that  14  MeV  neutrons  can  cause 
SEUs  in  large  memories  and  microprocessors,  the  impact  of  this  phenomena 
must  be  established  for  weapon  systems  that  have  to  be  hardened  to  radia¬ 
tion  from  thermonuclear  sources. 

The  event  rate  for  14  MeV  neutrons  ranges  from  one  upset  per  106  neu¬ 
trons/cm?  to  one  upset  per  5x10®  n/cm?  for  nMOS  dynamic  memories? 
and  between  one  upset  per  8x10*°  to  one  per  2x10*?  n/cm2  for  nM0S 
microprocessor.  Therefore,  SEUs  from  14  MeV  neutrons  can  occur  at 
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neutron  fluence  levels  well  below  the  levels  required  to  produce  signifi¬ 
cant  long  tern  displacement  damage. 

1-2  SCOPE. 

Three  basic  approaches  are  used  in  hardened  system  design: 

1)  Operation  through  a  Nuclear  Weapons  Environment  (NWE)  expo¬ 
sure  without  interruption  in  system  performance. 

2)  Operation  interrupted  but  restored  automatically  or  manually 
following  the  NWE  exposure. 

and 

3)  Functional  operation  maintained  through  the  NWE  exposure  with 
automatic  correction  for  radiation-induced  transient  effects 
(i.e.,  circumvention). 

In  the  first  case,  conventional  hardening  insures  that  no  critical  tran¬ 
sient  effects  will  occur  for  any  pulsed  ionizing  radiation  exposure  up  to  a 
critical  radiation  Intensity.  Systems  in  this  class  would  include  those 
hardened  for  operation  through  a  relatively  low-level  NWE  exposure,  or  a 
circumvented  system  for  exposures  up  to  the  system-defined  circumvention 
exposure  level.  The  critical  SEU  question  is  whether  it  is  possible  to  get 
a  sufficiently  high  neutron  exposure  for  the  pulsed  ionizing  radiation 
exposure  allowed  within  the  hardened  design. 

The  second  case  (i.e.,  automatic  or  manual  restoration  of  system 
operation)  has  not  been  considered  explicitly  in  this  study.  Neutron- 
Induced  SEUs  from  the  weapon  environment  may  increase  the  number  of  system 
Interruptions.  The  significance  of  the  Increase  will,  however,  depend  cri¬ 
tically  on  detailed  system  threat  scenarios  and  performance  requirements 
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not  considered  in  this  study.  Systems  of  specific  concern  could  include 
manned  aircraft  systems  operating  at  high  altitude. 


The  third  case  addresses  the  circumvented  system  for  which 
radiation-induced  errors  from  the  NWE  prompt  ionization  pulse  are  acknow¬ 
ledged  and  corrected.  The  critical  SEU  question  is  the  possibility  of 
neutron-induced  SEUs  occurring  after  the  system  assumes  normal  electrical 
operation  of  the  semiconductor  components.  The  concern  is  that  a  neutron 
pulse  of  sufficient  fluence  to  cause  SEUs  will  arrive  at  a  time  when  the 
system  assumes  normal  electrical  operation.  It  is  also  necessary  that  the 
overall  ionization  concomitant  with  the  neutron  pulse  is  insufficient  to 
re-initiate  system  circumvention. 

From  a  system  perspective  we  have  considered  systems  operated 
within  or  beyond  the  atmosphere.  If  both  the  system  and  weapon  are  endo- 
atmospheric,  the  scattering  of  neutrons  and  gamma  rays  are  important  as 
well  as  the  production  of  ionization  from  the  neutron  transport.  For  exo- 
atmospheric  systems  and  threats  we  only  have  to  consider  the  unscattered 
dispersion  of  the  NWE  prompt  ionization  and  neutron  pulses. 

1-3  SUMMARY  OF  RESULTS. 

The  susceptibility  of  circumvented  and  uncircumvented  endoatmos- 
pherlc  and  exoatmospherlc  systems  to  neutron  induced  single  event  upsets 
has  been  considered  in  very  general  terms  of  nuclear  weapon  threat  and  sys¬ 
tem  configuration. 

Fpr  both  exoatmospheric  and  endoatmospheric  uncircumvented  sys¬ 
tems  It  appears  that  transient  upset  effects  will  be  dominated  by  effects 
of  the  prompt  Ionization  pulse.  In  the  worst  case  neutron-induced  SEUs  are 
comparable  to  the  upset  levels  of  the  prompt  ionization  pulse.  For  current 
microelectronic  components  the  neutron-induced  SEUs  do  not  need  to  be 
explicitly  considered  for  uncircumvented  systems.  The  relative 


susceptibilities  should  be  reviewed  with  the  evolution  of  microcircuit 
technology,  particularly  for  possible  preferential  sensitivity  to  the 
localized  ionization  of  a  neutron-induced  reaction  compared  to  bulk  ioniza¬ 
tion  sensitivity. 

For  circumvented  systems  there  is  clearly  a  potential  considera¬ 
tion  for  exoatmospherlc  systems  and  a  possible  consideration  for  endoatmos- 
pherlc  systems.  Given  the  threat,  microcircuit  susceptibility  and  system 
performance  parameters  the  calculation  of  potential  concern  for  a  circum¬ 
vented  exoatmospheric  system  is  relatively  straightforward.  If  a  window  of 
susceptibility  exists,  two  immediate  hardening  approaches  can  be  con¬ 
sidered:  1)  increase  the  system  defined  electronics  recovery  time,  and/or 
2)  increase  the  sensitivity  of  the  system  radiation  circumvention  detec¬ 
tors.  The  acceptability  of  either  of  those  hardening  approaches  depends 
critically  on  many  system  performance  and  threat  parameters.  The  trade-off 
Issues  may  not  necessarily  be  easily  susceptible  to  resolution. 

The  high  atmospheric  concern  Is  comparable  to  the  exoatmospheric 
case.  In  some  cases  the  relative  neutron  transport  can  be  slightly 
enhanced  over  the  prompt  gamma  pulse  transport.  The  total  prompt  ioniza¬ 
tion  pulse  Is  also  reduced  by  the  increased  attenuation  of  the  low  energy 
photons.  The  relative  decrease  in  the  prompt  ionization  pulse  is,  however, 
relatively  unimportant  for  a  circumvented  system. 


SECTION  2 


UNCIRCUMVENTED  SYSTEM  -  EXOATMOSPHERIC 

When  both  the  system  and  weapon  are  exoatmospheric  (for  example, 
a  satellite  exposure  to  a  high  altitude  burst)  the  photon  and  neutron 
transport  are  without  attenuation.  The  prompt  gamma  ionization  dose  from  a 
"typical"  thermonuclear  weapon  Is  given  by, 

Y (prompt)  ■  rads(S1)/kiloton  (1) 

where  R  is  the  slant  range  between  the  weapon  and  the  system  in  kilo¬ 
meters.  The  prompt  dose  delivered  to  the  microcircuits  may  also  include  an 
additional  component  due  to  X-rays  which  depends  on  the  specific  weapon  and 
attenuation  of  the  external  X-ray  environment  by  the  system  structure.  In 
this  analysis  we  will  consider  only  the  prompt  gamma  component  which  will 
be  the  lower  limit  of  the  prompt  ionizing  dose.  Using  the  lower  limit  to 
the  prompt  ionizing  dose  will  give  an  upper  limit  to  the  expected  high 
energy  neutron  fluence  for  a  given  prompt  ionization  upset  level  of  the 
microcircuit. 

The  total  neutron  fluence  for  a  "typical"  thermonuclear  weapon  is 
given  approximately  as, 

*n  s  k6_y-0-l2  n/cmZ  -  kiloton  (E>0)  (2) 

where  R  Is  the  slant  range  in  kilometers.  In  this  analysis,  however,  we 
are  Interested  only  in  neutrons  of  sufficient  energy  to  cause  a  single 
event  upset  by  a  direct  n-a  reaction  with  a  silicon  atom.  Assuming  a  neu¬ 
tron  energy  threshold  of  4  MeV  for  the  silicon  n-a  reaction,  and  a  "typi¬ 
cal"  neutron  energy  spectrum  for  a  thermonuclear  weapon  as  shown  in  Table 
I,®  the  neutron  fluence  for  energies  greater  than  4  MeV  is  approximately 
0.167  of  the  total  fluence  so. 


Table  1.  Normalized  neutron  energy  outputs  for  different  energy  groups  for 
a  thermonuclear  burst. 8 


FRACTION  IN  GROUP  FOR 

GROUP  NEUTRON  ENERGIES  (MEV)  THERMONUCLEAR  WEAPON 


1 

1.6905+01 

- 

1.9640+01 

0.0 

2 

1.4918+01 

- 

1.6905+01 

0.0 

3 

1.4191+01 

- 

1.4918+01 

1.88714-02 

4 

1.3840+01 

- 

1.4191+01 

9.34254-03 

5 

1.2840+01 

- 

1.3840+01 

2.66169-02 

6 

1.2214+01 

- 

1.2840+01 

1.66622-02 

7 

1.1052+01 

- 

1.2214+01 

1.63678-02 

8 

1.0000+01 

- 

1.1052+01 

1.23974-02 

9 

9.0484+00 

- 

1.0000+01 

7.48258-03 

10 

8.1873+00 

- 

9.0484+00 

6.82320-03 

11 

7.4082+00 

- 

8.1873+00 

6.77521-03 

12 

6.3763+00 

- 

7.4082+00 

1.03201-02 

13 

4.9659+00 

- 

6.3763+00 

1.80706-02 

14 

4.7237+00 

• 

4.9659+00 

3.61700-03 

15 

4.0657+00 

- 

4.7237+00 

1.24302-02 

16 

3.0119+00 

- 

4.0657+00 

2.60380-02 

17 

2.3852+00 

- 

3.0119+00 

2.37305-02 

18 

2,3069+00 

- 

2.5852+OC 

3.74662-03 

19 

1.8268+00 

- 

2.3069+00 

2.56418-02 

20 

1.1080+00 

- 

1.8268+00 

6.44472-02 

21 

5.5023-01 

- 

1.1080+00 

8.84954-02 

22 

1.5764-01 

- 

5.5023-01 

9.13765-02 

23 

1.1109-01 

- 

1.5764-01 

1.16335-02 

24 

5.2475-02 

- 

1.1100-01 

1.10777-01 

25 

2.4788-02 

- 

5.2475-02 

5.40049-02 

26 

2.1875-02 

- 

2.4788-02 

5.68196-03 

27 

1.0333-02 

- 

2.1875-02 

9.26377-02 

28 

3.3546-03 

- 

1.0333-02 

1.16267-01 

29 

1.2341-03 

- 

3.3546-03 

7.38166-02 

30 

5.8295-04 

- 

1.2341-03 

2.32454-02 

31 

1.0130-04 

- 

5.8295-04 

2.02810-02 

32 

2.9023-05 

- 

1.0130-04 

1.90145-03 

33 

1.0677-05 

- 

2.9023-05 

0.0 

34 

3.0590-06 

- 

1.0677-05 

0.0 

35 

1.1254-06 

- 

3.0590-06 

0.0 

36 

4.1400-07 

- 

1.1234-06 

0.0 

37 

1.0000-11 

- 

4.1400-07 

0.0 

Total :  1.00000 
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(3) 


2.7  x  1011  _ / _ 2 

*n(>4  MeV)  s  n/cm  '  klloton 

The  expressions  for  the  prompt  ionizing  dose  and  high  energy  neutron 
fluence  are  plotted  in  Figure  1. 

Since  both  the  prompt  dose  and  high  energy  neutron  fluence  vary 
linearly  with  weapon  yield  and  are  both  inversely  proportional  to  R2,  their 
ratio  is, 

.  3  x  106  (4) 

Yp  rad(Si ) 

and  the  ratio  is  independent  of  weapon  yield  and  range.  Thus  if  the  total 
prompt  dose  was  1  rad(Si)  the  associated  high  energy  neutron  would  be  3xl06 
n/cm2  (E>4  MeV).  If  a  microcircuit  had  a  prompt  ionization  upset  level 
greater  than  1  rad(Si)  and  experienced  neutron-induced  single  event  upset 
at  3x1$  n/cm2  the  circuit  susceptibility  would  be  determined  by  the  neu¬ 
tron  SEUs. 


The  pulsed  ionization-induced  upset  level  in  microcircuits  is  a 
function  of  the  radiation  pulse  width.  For  some  microcircuits  the  circuit 
time  constants  are  long  compared  to  the  typical  radiation  pulse  width  and 
the  upset  level  can  be  defined  as  dependent  on  the  total  prompt  dose. 
Dose-dependent  microcircuit  types  include  dynamic  MOS  memories/micropro¬ 
cessors,  dynamic  bipolar  memories  and  CCD  arrays  for  radiation  pulse  widths 
of  less  than  ~1  ys.  Critical  total  dose  upset  levels  for  dynamic  MOS 
memories  are  on  the  order  of  0.1-10  rads(Si).  The  neutron  exposure  consis¬ 
tent  with  the  range  of  total  dose  upset  levels  in  dynamic  MOS  memories 
would  be  on  the  order  of  3  x  105  to  3  x  107  n/cm2  (E>4  Mev)  where  the 
higher  neutron  fluence  corresponds  to  the  least  sensitive  array  to  total 
dose-induced  upset. 
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Experimentally,  14  MeV  neutron-induced  upset  levels  measured  on 
dynamic  n-MOS  RAMs  vary  from  -108  -  5  x  109  n/cm2  for  4  k  memories, 

-107  -  5  x  109  n/cm2  for  16  k  memories,  and  -6  x  105  -  106  n/cm2  for 
early  64  k  memories. (1-5)  No  data  are  available  but  it  seems  reasonable 
to  assume  that  the  pulsed  ionization  total  dose  upset  on  the  large  memories 
is  also  relatively  low  (perhaps  0.1  -  1  rad(Si)). 

In  comparing  the  estimated  nuclear  weapon  high  energy  neutron 
fluence  and  prompt  ionization  dose  to  dynamic  RAM  data,  we  must  estimate 
the  relative  effects  of  the  4-14  MeV  weapon  spectrum  to  the  14  MeV  experi¬ 
mental  data  and  our  lower  bound  of  the  prompt  ionization  exposure.  The 
neutron  fluence  required  for  upset  increases  with  decreasing  energy  to  the 
4  MeV  threshold.  As  a  first-order  estimate  we  assumed  that  the  critical 
weapon  fluence  was  one-half  as  effective  as  the  14  MeV  monoenergetic 
exposure. 

The  comparison  of  neutron  and  prompt  ionization  upset  in  dynamic 
random-access  memories  is  summarized  in  Table  2.  For  the  4  k  and  16  k 
d-RAMs  it  appears  clear  that  the  upset  is  dominated  by  the  prompt  ioniza¬ 
tion  pulse.  The  prompt  ionization  and  neutron  SEU  levels  are  close  for  the 
64  k  RAM  but  additional  prompt  ionization  would  probably  clearly  push  the 
dominant  upset  to  that  of  the  prompt  ionization  pulse.  The  64  k  d-RAM  is  a 
sensitive  device  to  single  event  upset  and  it  might  be  expected  that  the 
SEU  rate  from  high  energy  protons  and  cosmic  rays  would  prohibit  applica¬ 
tion  In  a  satellite  system. 

The  pulsed  ionization  transient  upset  levels  for  static  semicon¬ 
ductor  memories  and  microcircuits  is  on  the  order  of  2  x  107  to  5  x  108 
rads  (S1)/s  for  pulse  widths  less  than  100  ns  and  varies  roughly  with  the 
square  root  of  the  pulse  width.  Assuming  weapon  pulse  width  of  20  ns,  the 
critical  prompt  gamma  dose  would  be  on  the  order  of  0.4  to  10  rads  (Si)  and 
the  associated  critical  neutron  upset  fluence  would  be  on  the  order  of  6  x 
103  -  3  x  105  n/cm2(>4  MeV).  There  is  very  little  data  on  the  neutron- 
induced  upset.  An  upset  level  of  -3  x  108  n/cm2 (E  *  14  MeV)  was  reported 


14  MeV  fluence 
for  SEU 
n/cm? 

Equiv.  NW 
fluence  >4  MeV 
n/cm2 

Estimated 

critical 

1  uc 

rads  (Si) 

$n (>4  Mev) 
Assoc,  with 

'  uc  , 
n/cnr 

4  k  d-RAM 

10?  -  5  x  109 

2  x  108  -  1010 

2 

6  x  106 

16  k  d-RAM 

107  -  5  x  109 

2  x  in7  -  io10 

0.5 

1.5  x  106 

64  k  d-RAM 

6  x  105  -  106 

1  -  2  x  106 

00.1 

3  x  105 

4  k  s-RAM 

-3  x  10? 

~6  x  108 

~2* 

6  x  106 

♦upset  level  is  pulse  width  dependent,  a  20  ns,  108  rads(Si)/s  pulse  was 
assumed. 


Table  2.  Comparison  of  RAM  measured  SEU  and  NWE  neutron  fluences. 
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for  the  Intel  4  k  (2147)  static  RAM.  The  comparison  of  prompt  ionization 
and  neutron  SEU  for  the  static  RAM  is  also  included  in  Table  2.  In  this 
case  it  appears  that  the  upset  is  clearly  dominated  by  the  prompt  ioniza¬ 
tion  pulse. 


In  summary,  it  appears  that  for  exoatmospheric  systems  and  NWE 
threats  that  microcircuits  will  upset  from  the  prompt  ionization  pulse  at 
exposure  levels  lower  than  those  required  for  neutron-induced  upset  for 
contemporary  digital  microelectronic  semiconductor  arrays.  Reported 
neutron-induced  upset  levels  for  advanced,  high  density  dynamic  semiconduc¬ 
tor  memories  are  near  upper  limit  estimated  NWE  levels  but  their  relative 
sensitivity  depends  on  the  trend  in  pulsed  ionization  upset  levels.  From 
this  analysis  it  appears  that  NWE  neutron-induced  single  event  upsets  do 
not  add  a  significant  susceptibility  to  digital  microelectronics  in  exoat¬ 
mospheric  systems  which  are  either  uncircumvented  or  circumvented  systems 
at  pulsed  ionizing  exposure  levels  below  the  circumvention  level. 
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SECTION  3 


CIRCUMVENTED  SYSTEM  -  EXOATMOSPHERIC 

The  next  case  of  concern  is  that  of  an  exoatmospheric  circum¬ 
vented  system  such  as  a  satellite  or  strategic  missile  in  flight.  We  will 
assume  that  the  prompt  ionization  pulse  is  of  sufficient  intensity  to  ini¬ 
tiate  system  circumvention.  For  the  purpose  of  this  discussion,  circumven¬ 
tion  can  be  described  as  a  system  design  such  that  the  radiation  pulse  is 
detected  and  transient  upsets  of  the  digital  microelectronics  are  allowed 
with  the  exception  of  critical  data  values.  After  a  predetermined  time,  it 
is  assumed  that  transient  effects  in  the  electronics  have  decayed  and  nor¬ 
mal  system  operation  can  be  resumed.  An  additional  time  interval  is 
required  to  update  the  system  and  complete  recovery.  The  circumvention 
time  Is  normally  defined  as  the  interval  between  initiation  and  recovery. 
The  critical  time  for  neutron-induced  single  event  upsets,  however,  is 
after  normal  electrical  performance  is  assumed  as  shown  in  Figure  2. 


initiation  electronics  system  recovery 

recovery 


Figure  2.  Assumed  timing  sequence  for  a  circumvented  system. 

Critical  parameters  for  neutron-induced  SEUs  in  a  circumvented 
system  are: 

1)  the  time  of  flight  of  the  neutron  pulse  compared  to  the  elec 
tronics  recovery  time  and, 

2)  the  high  energy  (>4  MeV)  neutron  fluence. 


The  neutron  time  of  flight  Is  shown  In  Figure  3  as  a  function  of  range  and 
neutron  energy.  Because  the  gamma  pulse  travels  at  the  speed  of  light,  it 
Is  reasonable  to  assume  that  the  delay  between  system  clrcumvehtlon  and  the 
leading  edge  of  the  neutron  pulse  Is  approximately  equal  to  the  14  MeV  neu¬ 
tron  time  of  flight. 

Using  the  1/R2  dispersion  for  neutrons  In  the  absence  of  scatter¬ 
ing  and  the  high  energy  neutron  yield  for  a  “typical"  thermonuclear  weapon 
(equation  3),  the  neutron  fluence  Is  shown  In  Figure  4  as  a  function  of 
tine  of  flight  (equivalent  to  range)  for  weapons  of  different  yield. 

The  Influence  of  possible  system  parameters  Is  Illustrated  In 
Figure  5.  If  we  assume  that  the  electronics  recovery  time  Is  5  ms,  then 
neutron-induced  SEUs  occurring  earlier  will  be  corrected  with  the  rest  of 
the  scrambled  electronics  memory.  If  we  also  assume  that  the  microcircuit 
sensitivity  is  3  x  1$  n/cm2  (E>4  MeV),  which  is  on  the  order  of  that 
observed  for  an  n-MOS  static  RAM,  neutron  exposures  below  the  critical 
upset  level  will  not  be  of  concern  even  if  the  neutron  pulse  arrives  later 
than  the  electronics  recovery  time.  The  window  of  vulnerability  for  the 
system  to  neutron-induced  SEUs,  for  the  assumed  parameters  would  be  for 
thermonuclear  weapon  exposures  of  yield  greater  than  100  kT  at  ranges 
greater  than  approximately  200  km.  We  feel  that  a  window  of  vulnerability 
can  exist  for  circumvented  exoatmospheric  systems.  The  degree  of 
existence,  or  absence,  of  the  window  of  vulnerability  depends  strongly  on 
the  specific  system  performance  characteristics  and  threat  and  no  general 
conclusions  can  be  made. 


Unscattered  Neutron  Traversal  Time,  seconds 


Range,  km 


Figure  3.  Unscattered  neutron  traversal  time  for  various  monoenergetic 
source  energies  as  a  function  of  range. 
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Figure  5.  Possible  exoatmospheric  window  of  susceptibility. 
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SECTION  4 


ENDOATMOSPHERIC  SYSTEMS 

In  the  analysis  for  the  exoatmospheric  environment,  the  high 
energy  neutron  SEU  phenomena  was  addressed  by  merely  calculating  the  spher¬ 
ical  divergence  of  the  prompt  gamma  rays  and  the  neutrons.  Two  cases  of 
Interest  were  established,  a)  the  system  does  not  circumvent  and  a  suffi¬ 
cient  number  of  high  energy  neutrons  arrive  at  the  target  to  cause  SEUs. 

In  this  case  the  arrival  time  of  the  neutrons  is  of  no  interest.  And  b) 
the  system  circumvents  but  a  sufficient  number  of  high  energy  neutrons  to 
cause  SEUs  arrive  at  the  target  after  it  has  recovered.  In  this  case  the 
neutron  arrival  times  must  be  known.  If  both  of  these  cases  are  considered 
for  the  endoatmospheric  environment,  then  the  attenuation  of  the  prompt 
gamma  rays,  the  attenuation  and  energy  degradation  of  the  neutrons  and  the 
time  dispersion  of  the  neutrons  must  be  calculated  for  transport  through 
non-uniform  air. 

Neutron  and  gamma  ray  transport  through  air  has  been  studied 
extensively  and  a  large  volume  of  calculated  data  is  avai lable8-12.  The 
air  transport  problem  has  been  addressed  with  a  one  dimensional  transport 
code,  ANISN,  in  uniform  air  and  with  Monte  Carlo  calculations  in  both  uni¬ 
form  and  non-uniform  air.  Comparisons  between  calculations  made  in  uniform 
and  nonuniform  air  Indicate  that  over  a  wide  range  of  altitudes  and  slant 
ranges,  the  uniform  air  calculations  can  be  applied  to  any  endoatmospheric 
scenario  by  using  the  concept  of  mass  Integral  scaling  (MIS).  This  concept 
means  that  the  transport  of  gammas  and  neutrons  is  simply  a  function  of  how 
much  air  mass  the  radiation  must  pass  through  from  source  to  target.  In 
addition  to  the  attenuation  through  a  given  air  mass,  spherical  divergence 
(1/4* R2)  must  be  taken  Into  account  In  order  to  determine  the  prompt 
gamma  dose  rate  and  neutron  fluence  for  a  given  range  and  air  mass.  Thus 
the  radiation  environment  reaching  the  target  is  determined  from  a  relation 
of  the  form 


where  ENVy  is  the  environment  of  interest  at  the  target,  S  is  the  source 
strength,  R  the  range,  y  an  attenuation  coefficient  in  cm^/g  and  ma  the 
air  mass  in  g/cm?. 

If  the  results  of  calculations  made  in  uniform  air  are  used  to 
determine  the  value  of  y  and  MIS  is  applied,  the  results  should  be  within  a 
factor  of  2  for  most  scenarios.  Thus,  to  determine  the  extent  of  the  high 
energy  neutron  SEU  phenomena  in  air,  uniform  air  calculations  will  be  used. 
The  main  reason  for  using  the  uniform  air  results  is  that  the  available 
literature  for  uniform  air  calculations  includes  neutron  energy  and  time 
dependence  whereas  the  nonuniform  air  results  include  only  the  integrated 
neutron  results  over  the  whole  energy  range.  For  the  SEU  phenomena,  only 
the  neutrons  with  energies  above  -4  MeV  are  of  interest. 

The  endoatmospherlc  scenario  can  be  extrapolated  from  the 
exoatmospheric  results  by  including  an  exponential  attenuation  factor  for  a 
given  air  mass.  Thus  for  a  given  scenario  In  air  the  range  and  total  air 
mass  must  first  be  determined. 

4-1  CALCULATION  OF  AIR  MASS. 

The  density  of  air  as  a  function  of  altitude  is  assumed  to  follow 
an  exponential  decrease  from  the  sea  level  value  of  1.22  mg/cm3.  Using 
this  relation  graphs  and  nomographs  have  been  established  to  determine  the 
air  mass  for  a  given  range  and  altitude  of  source  and  target.  For  the  case 
where  the  source  is  exoatmospheric  and  directly  above  the  target,  the  air 
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mass  between  source  and  target  can  be  determined  for  a  given  target  alti¬ 
tude  from  Figure  6.13  in  this  case  a  simple  integration  of  the  variable 
density  from  the  target  altitude  to  »  has  been  performed.  For  the  case 
where  the  source  and  target  are  at  the  same  altitude,  the  air  mass  can  be 
determined  from  Figure  7.13  jn  this  figure  the  coaltitude  is  given  on 
the  vertical  scale  and  the  curves  for  coaltitude  range  are  given  for  .5  km 
to  40  km. 


If  the  source  and  target  are  not  at  the  same  altitude,  then  the 
nomograph  given  in  Figure  8  must  be  used. 13  Directions  for  using  the 
nomograph  are  given  in  the  figure. 

Once  the  range  and  air  mass  are  determined  for  a  given  scenario, 
then  the  prompt  gamma  dose  rate  and  high  energy  neutron  fluence  can  be 
determined  from  the  spherical  divergence  and  uniform  air  transport  calcula¬ 
tions. 


4-2  CALCULATION  OF  PEAK  DOSE  RATE  IN  AIR. 

The  maximum  dose  rate,  which  can  result  in  ci rcumventon,  will  be 
due  to  the  prompt  gamma  rays.  Attenuation  of  the  prompt  gammas  in  air,  due 
to  the  interaction  of  the  gammas  with  air  molecules  is  given  by  a  total 
mass  attenuation  coefficient  and  a  total  mass  absorption  coefficient.  Both 
of  these  exponential  decay  factors  are  energy  dependent.  The  mass  attenua¬ 
tion  coefficient  assumes  that  all  gammas  which  undergo  an  interaction  in 
air  are  removed  from  the  beam.  Hence,  coherent  scattering  is  not  accounted 
for.  The  mass  absorption  coefficient  is  determined  by  including  coherent 
scattering  and  hence  the  removal  cross  section  is  somewhat  lower.  Table  3 
is  a  list  of  the  mass  absorption  and  mass  attenuation  coefficients  for 
gammas  of  various  energies  in  dry  air. 
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Figure  8.  Computation  of  mass  penetrated  between  two  points  in  the  atmosphere 
1  kg/m2  =  .1  g/cm2.13 


Table  3.  Total  mass-absorption  and  mass-attenuation  cross  sections 
in  air  for  various  gamma-ray  energies. 


Gamma-Ray 

Energy 

(MeV) 

— 

Total  Mass-Absorption 

Cross  Section 
(cm2  /g) 

Total  Mass-Attenuation 
Cross  Section 
(cm2 /g) 

0.5 

2. 95{ -2) 

8.65( -2) 

1.0 

2.82-2 

6.26(-2) 

2.0 

2. 37( -2) 

4. 38( -2) 

4.0 

1.92(-2) 

4. 06( -2) 

6.0 

1. 72( —2) 

2.49(-2) 

8.0 

1. 60( -2) 

2. 1 9(  -  2) 

10.0 

1. 54( -2) 

2. 0 1 ( -2) 

Note:  Mass  fractions  =  0.755  nitrogen,  0. 

232  oxygen,  0.013  argon. 

For  the  purposes  the  present  calculations,  the  mass  absorption 
coefficient  will  be  used  since  it  represents  a  more  realistic  case.  A  plot 
of  the  attenuation  of  the  prompt  gammas  as  a  function  of  air  mass  is  given 
in  Figure  9  for  0.5,  2  and  10  MeV  gammas. 

The  dose  rate  from  prompt  gammas  at  the  target  can  be  determined 
if  the  following  information  is  known. 

a)  Air  mass 

b)  Range 

c)  Source  strength  and  energy  spectrum  for  weapons  type  -  gammas 
per  kiloton  and  energy  distribution. 

d)  Emission  time  of  gammas 


To  convert  the  gamma  flux  to  dose  rate  in  silicon,  the  Henderson  flux  to 
dose  conversion  factor  is  used.  For  the  exoatmospheric  calculations  the 


dose  per  klloton  versus  range  was  determined  from  Figure  1.  This  figure 
was  derived  assuming  a  "typical"  thermonuclear  source  and  using  the 
Henderson  gamma  to  dose  (Si)  conversion  for  1  MeV  gammas. 13 


4-3  HIGH  ENERGY  NEUTRON  TRANSPORT  IN  AIR. 

For  the  high  energy  neutron  SEU  phenomena,  the  high  energy  neu¬ 
tron  fluence  must  be  determined  at  the  target  after  transporting  through  a 
given  air  mass.  Therefore,  not  only  the  total  neutron  attenuation  must  be 
known,  but  the  degradation  of  the  energy  spectrum  as  well.  St  raker  has 
calculated  the  angular  and  energy  dependence  of  neutron  transport  in  uni¬ 
form  dry  air  from  both  discrete  neutron  sources  and  a  typical  thermonuclear 
source. 8  The  one-dimensional  transport  code  ANISN  was  used.  Calcula¬ 
tions  were  carried  out  to  4800  meters  (533  g/cm2)  for  the  12.2-15  MeV 
source  and  out  to  1800  meters  (200  g/cm2)  for  the  thermonuclear  source. 

Using  Straker's  results,  the  4nR2  neutrons  per  source  neutron 
values  were  calculated  for  neutron  energies  greater  than  4.06  MeV  at  the 
target  as  a  function  of  air  mass.  The  results  are  shown  in  Figure  10  for 
the  12.2-15  MeV  source  anc.  che  thermonuclear  source.  The  energy  spectrum 
of  the  thermonuclear  source  is  shown  in  Table  4. 

Also  shown  In  Figure  10  is  the  attenuation  of  1  MeV  prompt  gammas 
using  the  mass  absorption  coefficient.  The  results  of  the  neutron  trans¬ 
port  calculations  indicate  that  for  air  masses  of  greater  than  100-200 
g/cm2  the  attenuation  of  neutrons  can  be  represented  fairly  well  by  a 
constant  exponential  attenuation  factor  yn  in  g/cm2.  From  the  calcu¬ 
lated  results  in  Figure  10,  neutron  attenuation  coefficients  (removal  cross 
sections)  were  determined.  These  are  presented  in  Table  5  for  the  12.2-15 
MeV  and  thermonuclear  sources. 
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Table  4.  Neutron  energy  spectrum  for  a  thermonuclear  source. 


Energy  (MeV) 

Fraction  in  Group 

12.2-15 

.0706 

10.-12.2 

.0256 

8.19-10. 

.0141 

6.36-8.19 

.0147 

4.97-6.36 

.0180 

4.07-4.97 

.0170 

3.01-4.07 

.0260 

2.46-3.01 

.0190 

2.35-2.46 

.0050 

1.83-2.35 

.0280 

1.11-1.83 

.0620 

.55-1.11 

.0850 

.55  < 

.615 

Table  5.  Neutron  attenuation  coefficients  versus  air  mass  for  12.2-15  MeV 
and  thermonuclear  sources. 


12.2-15  MeV 

Source 

Thermonuclear 

Source 

(En  2  4-07  MeV 

is  16%) 

Air  mass 

cm^/g 

Air  mass 

jfn 

c-r/g 

g/cm? 

g/cm2 

33.3 

.0145 

33.3 

.070 

66.6 

.0221 

66.6 

.0495 

133 

.0283 

100 

.0434 

200 

.0311 

133 

.0418 

266 

.0329 

200 

.0400 

400 

.0348 

533 

.0361 
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Although  the  value  of  pn  is  not  constant  over  the  entire  range 
of  air  masses  of  interest  for  most  scenarios,  it  is  fairly  constant  above 
100  g/cm2.  A  nominal  value  of  .033  cm2/g  could  be  used  for  the  12.2-15 
MeV  source  and  .040  cm2/g  for  the  thermonuclear  source. 

4-4  TIME  DISPERSION  OF  NEUTRONS. 

For  the  SEU  case  where  the  system  circumvents,  the  neutron  arri¬ 
val  times  must  be  known.  For  an  exoatmospheric  environment,  the  neutron 
arrival  time  is  determined  from  the  direct,  unscattered  neutron  time  of 
flight.  However  for  transport  in  air  only  a  few  of  the  neutrons  reach  the 
target  without  undergoing  either  an  elastic  or  inelastic  collision.  There¬ 
fore,  the  arrival  becomes  dispersed  in  time.  StrakerlO  has  performed 
time-dependent  Monte-Carlo  transport  calculations  in  uniform  air  for  a 
source  50  ft.  above  ground.  Results  are  tabulated  for  the  4uR2  neutron 
flux  as  a  function  of  energy  and  time  for  slant  ranges  up  to  1300  meters. 
Neutron  sources  were  discrete  energy  bins  and  a  fission  spectrum.  No  cal¬ 
culations  were  made  for  a  thermonuclear  source.  In  order  to  determine  the 
time  dispersion  for  the  high  energy  neutrons  of  interest  to  SEUs,  calcula¬ 
tions  of  the  total  4irR20(En  2  4*06  MeV)  have  been  made  versus  time  for 
12.2-15  MeV  source.  The  results  are  shown  in  Figure  11  for  slant  ranges  of 
700  and  1300  meters.  These  ranges  correspond  to  77  and  143  g/cm2  air 
mass.  The  two  histogram  plots  in  Figure  11  show  the  fraction  of  the  total 
high  energy  neutron  fluence  reaching  the  detector  for  the  various  time 
bins.  Also  shown  (shaded  area)  Is  the  time  bin  bounded  by  the  direct  tran¬ 
sit  time  for  the  15  MeV  neutrons  to  the  near  side  of  the  detector  and  the 
transit  time  for  12.2  MeV  neutrons  to  the  far  side  of  the  detector.  For 
the  700  meter  slant  range,  the  detector  extends  from  675  to  725  meters  and 
for  the  1300  meter  slant  range,  the  detector  extends  from  1225  to  1375 
meters.  Thus  the  shaded  area  represents  the  maximum  range  of  time  of 
flight  for  the  unperturbed  source  neutrons.  If  the  nominal  time  of  flight 
for  the  time  dispersed  neutrons  is  compared  to  the  time  of  flight  for 
unperturbed  neutrons  then  the  arrival  times  are  seen  to  approximately 
double  for  both  slant  ranges. 
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Figure  11.  Time  dispersion  of  high  energy  neutrons 


SECTION  5 


COMPARISON  OF  EXOATMOSPHERIC  AND  ENDOATMOSPHERIC  ENVIRONMENTS 

In  the  section  on  the  exoatmospheric  environment,  two  sets  of 
curves  were  developed  to  determine  the  high  energy  neutron  fluence  from  a 
typical  thermonuclear  weapon  for  the  two  cases,  a)  no  circumvention,  and  b) 
circumvention  with  neutrons  arriving  after  the  system  has  recovered.  Using 
the  concept  of  exponential  attenuation  in  addition  to  spherical  divergence, 
the  exoatmospheric  results  can  be  extended  to  the  endoatmospheric  environ¬ 
ment  If  the  amount  of  air  mass  between  source  and  target  is  known. 

For  the  case  where  the  system  does  not  circumvent,  i.e.,  the  peak 
dose  rate  from  prompt  gammas  is  less  than  the  threshold  level  for  circum¬ 
vention,  the  fluence  of  high  energy  neutrons  reaching  the  target  is  reduced 
from  the  exoatmospheric  calculation  by  the  ratio  between  high  energy  neu¬ 
tron  attenuation  and  the  prompt  gamma  ray  attenuation. 

For  example,  if  the  air  mass  is  300  g/cm2,  corresponding  to  the 
target  at  26,000  ft.  and  the  source  is  exoatmospheric  and  directly  above 
the  target,  then  the  ratio  between  the  thermonuclear  high  energy  neutron 
attenuation,  exp(-vnma),  and  the  1  MeV  prompt  gamma  attenuation, 
exp(-viYma),  is  .029.  Thus  the  fluence  of  high  energy  neutrons  in  300 
g/cm2  of  air  would  be  only  .029  times  the  exoatmospheric  fluence  for  the 
same  circumvent  threshold  level. 

Closer  inspection  of  Figure  10  and  the  data  in  Table  5  indicate 
that  this  simplified  procedure  cannot  be  applied  to  air  masses  of  less  than 
about  150  g/cm2.  For  air  masses  less  than  this  value  the  high  energy 
neutron  attenuation  coefficient  Is  not  constant.  In  fact  the  results  for 
the  12.2-15  MeV  source  attenuation  curve  lies  slightly  above  the  1  MeV 
gamma  attenuation  curve  below  100  g/cm2  air.  The  thermonuclear  curve  is 
displaced  downward  since  only  16%  of  the  source  neutrons  have  energies  >  4 
MeV  to  begin  with. 
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Based  on  these  results,  the  following  conclusions  can  be  drawn 
for  the  endoatmospheric  environment  for  the  case  where  the  system  does  not 
ci rcumvent : 

1)  For  air  masses  up  to  -150  g/cm2  the  high  energy  neutron 
fluence  will  be  essentially  the  same  as  for  the  exoatmospher- 
ic  environment.  Thus,  Figure  1  can  be  used  for  the  typical 
thermonuclear  source. 

2)  For  air  masses  greater  than  150  g/cm2,  the  high  energy  neu¬ 
tron  fluence  will  be  lower  than  for  the  exoatmospheric 
environment.  The  reduction  factor  can  be  estimated  from  the 
ratio 


The  second  case  to  consider  is  the  one  for  which  the  source  is  far  enough 
away  such  that  the  high  energy  neutrons  arrive  after  the  system  has 
recovered  from  a  circumvention.  For  this  case  the  time  dispersion  of  the 
neutrons  in  the  atmosphere  as  well  as  the  attenuation  must  be  determined. 
The  arrival  time  for  unperturbed  neutrons  versus  range  is  given  in  Figure 
3  for  various  monoenergetic  neutrons. 

For  a  system  off  time  of  only  1  ms  the  range  for  4  MeV  to  14  MeV 
unscattered  neutrons  is  28  to  52  km.  The  results  of  the  time  dispersion 
calculations  given  In  Figure  9  indicate  that  only  30-40%  of  the  high  energy 
neutrons  arrive  at  times  greater  than  twice  the  unscattered  arrival  times 
for  ranges  up  to  1.3  km.  Unfortunately,  no  tablulated  data  in  the  open 
literature  could  be  found  for  calculations  of  time  dispersion  at  greater 
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ranges.  However  it  can  be  assumed  that  the  time  dispersion  is  more  a  func¬ 
tion  of  air  mass  than  range.  The  results  of  Figure  11  extend  to  143 
g/cm2  air  mass.  Even  If  a  factor  of  two  increase  in  dispersion  is 
assumed  with  only  a  factor  of  3  decrease  in  high  energy  neutron  fluence 
then  the  range  can  only  be  reduced  to  14-26  km  for  the  high  energy  neutrons 
for  a  1  ms  arrival  time.  For  a  14-26  km  range  the  minimum  air  mass  for  an 
endoatmospheric  scenario  would  be  with  the  source  exoatmospheric  and  the 
target  at  some  altitude  directly  below  the  source.  If  we  assume  an  alti¬ 
tude  of  50,000  ft.  for  the  target,  the  air  mass  would  be  about  120  g/cm2 
from  Figure  6.  Using  the  attenuation  calculations  ^  Figure  10  for  the 
12.2-15  MeV  neutron  component,  the  attenuation  of  12u  g/cm2  is  .038  or  a 
factor  of  26  lower  than  the  corresponding  exoatmospheric  calculation.  For 
a  coaltitude  scenario  at  50,000  ft.  the  air  mass  for  the  coaltitude  range 
of  14-26  km  is  230-450  g/cm2.  This  would  correspond  to  a  6x10-4  to 
10“7  reduction  in  high  energy  neutrons  from  the  exoatmospheric  calcula¬ 
tion. 

Based  on  the  endoatmospheric  calculations  for  the  case  where  the 
system  does  circumvent  the  following  conclusions  can  be  drawn. 

1)  The  high  energy  neutron  fluence  will  be  less  than  the  high 
energy  neutron  fluence  for  the  exoatmospheric  environment. 

2)  The  time  dispersion  in  the  atmosphere  will  not  be  great 
enough  to  significantly  reduce  the  required  range  between 
source  and  target. 

3)  The  minimum  amount  of  air  mass  for  the  required  range  will 
occur  for  the  scenario  where  the  source  is  exoatmospheric  and 
directly  above  the  high  altitude  target. 

4)  If  both  source  and  target  are  located  in  the  atmosphere,  then 
the  air  mass  for  the  required  range  will  significantly  reduce 
the  high  energy  neutron  fluence. 


SECTION  6 
CONCLUSIONS 

Concern  for  neutron-induced  single  event  upsets  compromising  the 
performance  of  hardened  military  systems  can  be  listed  in  order  of  decreas¬ 
ing  concern  as  follows: 

1)  Circumvented,  exoatmospheric  systems. 

2)  Circumvented,  high  endoatmospheric  systems. 

3)  Uncircumvented,  high  endoatmospheric  systems. 

4)  Uncircumvented,  exoatmospheric  systems. 

5)  Low  endoatmospheric  systems. 

The  degree  of  concern  will  be  a  function  of  detailed  system  performance 
requirements,  detailed  threat  scenarios,  and  microcircuit  SEU  susceptibi¬ 
lity.  Subjectively,  the  concern  for  circumvented,  exoatmospheric  systems 
seems  clear  for  systems  using  modern  MOS  or  bipolar  semiconductor  memories 
(e.g.,  >  1  k  bit/chip).  The  concern  for  circumvented,  high  endoatmospheric 
systems  Is  not  as  clear  but  approaches  that  of  exoatmospheric  systems.  The 
concern  for  uncircumvented  and  low  endoatmospheric  systems  is  not  clear  and 
probably  represents  no  significant  problem  with  today's  semiconductor  com¬ 
ponents. 

Parameters  of  circumvented  system  can  be  adjusted  to  minimize  or 
eliminate  a  neutron  SEU  problem  with  trade-offs  in  overall  system  perfor¬ 
mance.  For  example,  the  defined  system  recovery  time  can  be  increased 
and/or  the  radiation  detector/circumvention  level  can  be  decreased  (i.e., 
Increased  in  sensitivity).  If  the  recovery  time  is  increased,  the  total 
circumvention  time  will  increase  and  the  increased  "off"  time  may  degrade 
overall  system  performance.  If  the  circumvention  level  is  decreased,  the 
number  of  circumventions  for  complex  threat  scenarios  may  increase,  also 
Increasing  the  system  "off"  time  with  possible  degradation  in  system  per¬ 
formance. 
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An  alternative  approach  to  system  hardening  is  to  employ  micro- 
circuits  of  acceptable  SEU  susceptibility.  Unfortunately,  the  trend  seems 
to  be  an  increase  in  SEU  susceptibility  as  we  go  toward  more  complex 
advanced  semiconductor  devices.  For  the  specific  concern  of  weapon-envi¬ 
ronment  effects  we  can  at  least  be  confident  that  the  microcircuit  suscep¬ 
tibility  cannot  increase  above  the  limit  defined  by  errors  induced  by  cos¬ 
mic  rays  and  high  energy  protons.  The  cosmic  ray  and  high  energy  proton 
SEU  susceptibility  is  a  basic  operational  concern  for  satellite  systems  and 
cosmic  ray  susceptibility  is  a  potential  concern  for  high  altitude 
avionics. 

On  balance,  neutron-induced  single  event  upsets  represent  an 
additional  hardening  concern  for  some  classes  of  military  systems  but  do 
not  represent  a  major  concern  for  the  broad  range  of  military  systems  and 
does  not  represent  a  significant  exploitable  vulnerability  to  existing 
systems. 
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APPENDIX 

LUMPED  MOOEL  ANALYSIS  OF  TWO-DIMENSIONAL  CARRIER  COLLECTION 


It  is  generally  assumed  that  the  transient  current  of  a  reverse- 
biased  p-n  junction  from  a  single  particle  ionization  track  is  a  narrow 
pulse.  If  the  pulse  width  is  short  compared  to  the  smallest  circuit  time 
constant  then  the  circuit  response  will  be  determined  by  the  total  collected 
charge  and  the  details  of  the  waveform  are  not  important.  On  the  other 
hand,  if  the  pulse  width  is  comparable  to  the  circuit  time  constants  details 
of  the  waveform  will  be  important  in  establishing  the  circuit  upset  level. 

Charge  generated  by  the  particle  ionization  track  in,  and  near, 
the  junction  depletion  layer  will  be  collected  quickly.  Detailed  one¬ 
dimensional  calculations  for  realistic  doping  profiles  and  including  effects 
of  hole  and  electron  transport  (i.e.,  funnel ing)  predict  a  subnanosecond 
current  pulse  (Ref.  1-3).  This  fast  current  pulse  is  dominated  by  carrier 
collection  in  the  depletion  layer  which  is  enhanced  by  funnel ing. 

Carriers  generated  deeper  in  the  bulk  semiconductor  will  diffuse 
to  the  junction  with  some  time  delay  before  collection.  If  the  bulk  semi¬ 
conductor  were  ionized  uniformly  this  would  be  the  familiar  diffusion  compo¬ 
nent  of  photocurrent.  The  diffusion  photocurrent  decays  with  a  time  con¬ 
stant  on  the  order  of  the  minority  carrier  lifetime.  For  modem  non-gold- 
doped  microcircuits  the  high  resistivity  substrate  lifetime  can  be  on  the 
order  of  microseconds.  Time  constants  on  this  order  would  be  significant  In 
the  response  of  virtually  all  types  of  microcircuits. 

Intuitively,  the  diffusion  current  resulting  from  a  single  parti¬ 
cle  Ionization  will  be  less  than  that  observed  by  uniform  carrier  generation 
because  of  carrier  diffusion  away  from  the  Ionization  track  (Figure  1). 

This  would  suggest  a  photocurrent  pulse  which  is  smaller  and  narrower  than 
the  uni  form- ionization  photocurrent  (Figure  2). 
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Figure  1.  Lateral  diffusion  of  single  particle-generated  carriers. 
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Figure  2.  Possible  photocurrent  waveforms. 
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The  purpose  of  this  analysis  is  to  estimate  the  extent  of  the 
lateral  carrier  diffusion  with  the  use  of  a  simple  two-dimensional  model  of 
carrier  flow  and  collection. 

A  two-dimensional,  cyl Indr ically- symmetric  bulk  semiconductor  re¬ 
gion  Is  shown  In  Figure  3.  For  this  analysis  we  will  approximate  carrier 
densities  by  their  average  value  over  finite  volumes.  If  a  large  number  of 
small  volumes  were  assumed  this  would  simply  be  a  finite-difference  repre¬ 
sentation  of  the  continuous  equations  of  carrier  distribution  and  flow.  In 
this  case,  however,  we  know  that  the  carrier  distribution  will  vary  strongly 
with  distance  near  the  junction  and  Ionization  track  and  more  slowly  with 
distance  far  into  undisturbed  bulk  semiconductor.  We  therefore  can  select 
smaller  volumes  near  the  junction  and  ionization  track  and  Increase  the  vol¬ 
umes  as  we  move  out.  In  this  way  we  can  maintain  reasonable  accuracy  In  the 
analysis  with  a  much  simpler  model. 

For  the  region  to  be  analyzed,  as  shown  In  Figure  3,  we  will  as¬ 
sume  that  the  reverse-biased  p-n  junction  in  In  the  x-r  plane  at  x  =  0  and 
the  particle  Ionization  track  Is  down  the  center  at  r  -  0. 

The  critical  first-order  semiconductor  equations  to  be  differenced 
are  those  relating  minority  carrier  density,  generation,  recombination  and 
diffusion  In  two  dimensions.  Assuming  uniformly  doped  n-type  material  for 
the  lumped  volume  at  x  *  1  and  r  *  j  we  have, 

,«co«b.1ij  ‘  q,1.J  <pi.j  -  "n’/'p  ,U 

for  the  carrier  recombination  where  j  Is  the  lumped  element  volume,  p^  ^ 
Is  the  average  minority  carrier  density,  p^  Is  the  thermal  equilibrium 
minority  carrier  density,  Is  the  minority  carrier  lifetime.  If  we  assume 
that  the  particle- Induced  carrier  density  Is  large  compared  to  the  thermal 
equilibrium  value, 
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Carrier*  generation  In  the  lumped  volume  can  be  expessed  as, 
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where  go(t)  Is  the  Induced  time-dependent  carrier  generation  rate.  In  this 
analysis,  this  term  will  be  Included  only  in  the  lumped  elements  along  the 
x-axis. 
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Increase  and  decrease  in  carrier  density  in  the  lumped  element  can 
be  represented  by  a  carrier  storage  term  as. 


*  storage,  j  *  (<l  V1 ,  j  ’  ‘  dpi .  j/dt> 


(4) 


Carrier  diffusion  may  be  either  axial  or  radial.  For  this  term 
the  current  flow  Is  determined  by  the  difference  in  average  value  of  carrier 
density  in  contiguous  elements  and  the  spacing  between  the  central  points  of 
the  elements.  Axial  diffusion  current  is  then  approximately: 
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between  the  1th  and  ( 1+1 )th  elements  where  Ax  is  the  area  of  the  axial 
Interface  and  Op  Is  the  minority  carrier  diffusion  constant.  Similarly, 
for  the  radial  current. 
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where  V  is  the  applied  junction  voltage  and  the  boundary  condition  is  con¬ 
sidered  uniform  across  the  junction  plane. 

Reviewing  Equation  1-6  each  involves  constants  of  the  bulk  semi¬ 
conductor,  area,  and  lumped  geometry  multiplied  by  terms  of  the  minority 
carrier  density.  Linvill  defined  these  collective  constants  be  terms  and 
symbols  (Ref.  4).  Thus  a  lumped  model  network  could  be  defined  in  terms  of 
carrier  density,  current  and  voltage  as  shown  in  Figure  4  for  a  one-dimen¬ 
sional  model.  In  general,  the  lumped-model  can  be  used  to  represent  carrier 
flow  by  drift  as  well  as  diffusion. 

The  lumped  model  network  generally  has  fewer  nodes  than  a  finite- 
difference  model  but  often  still  has  too  many  for  convenient  closed  form 
solution.  A  wide  variety  of  circuit  analysis  computer  programs  is  available 
for  the  analysis  of  linear  and  non-linear  networks.  Unfortunately,  the 
electrical  circuit  network  codes  deal  only  in  terms  of  current  and  voltage 
while  the  general  lumped  model  network  must  Involve  current,  voltage  and 
carrier  density  [NET-2  has  a  provision  of  direct  analysis  of  a  limited  set 
of  lumped  model  networks  but  this  provision  Is  not  widely  used  or  easily 
available]  (Ref.  5). 

For  this  analysis,  the  difficulty  can  be  resolved  because  we  will 
assume  that  the  junction  Is  always  reverse  biased.  The  excess  minority  car¬ 
rier  density  Is  then  fixed  at  zero  at  the  junction  and  the  actual  junction 
voltage  Is  not  a  necessary  parameter  for  the  analysis.  Now  the  only  param¬ 
eters  of  the  lumped  model  analysis  are  minority  carrier  density  and  cur¬ 
rent.  An  electrical  network  can  now  be  defined  by  considering  the  minority 


Figure  4.  One-dimensional  lumped  model. 


carrier  density  as  a  “voltage".  The  junction  photocurrent  will  be  the  cal¬ 
culated  current  that  flows  through  a  short-circuit  at  the  junction.  An 
electrically  equivalent  network  for  a  two-dimensional  model  Is  shown  In  Fig¬ 
ure  5. 


The  lumped  model  of  Figure  5  approximates  the  bulk  semiconductor 
volume  by  five  regions  along  the  x-axIs  and  five  radial  volumes.  The  geom¬ 
etry  was  scaled  in  units  of  the  minority  carrier  diffusion  length  Lp.  The 
points  of  average  carrier  density  for  the  axial  and  radial  dimensions  were 
selected  as  shown  In  Table  I.  As  discussed  previously  the  points  are  spaced 

Table  I.  Coordinates  for  Average  Values  of  Carrier  Density. 


Xi  s  O*1  LP 

X2  =  0.2  Lp 

x3  *  0.3  L 
3  P 

x,.  =  0.5  L 
^  P 

Xc  =  0.8  L 

^  P 


rl  =  0*1  Lp 

r2  *  0.2  L 
L  P 

r»  =  0.3  L 
3  P 

r4  =  0.5  L 
H  P 

r5  =  0.8  L 
3  P 


closer  near  the  junction  and  center  axis  to  preserve  accuracy  with  a  rela¬ 
tively  small  number  of  lumped  regions  (Ref.  6).  The  volumes  associated 
with  each  lumped  region  are  presented  in  Table  II.  In  each  as  we  have. 
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In  terms  of  the  circuit  network  the  value  of  the  carrier  recombination  "con¬ 
ductances"  Is 


Figure  5.  Two  dimensional  equivalent  lumped  model. 


Table  II.  Lumped  Model  Volumes  in  Units  of  nL3 


r  0 

1 

2 

3 

4 

X 

0.001 

0.003 

0.005 

0.016 

0.039 

1 

0.001 

0.003 

0.005 

0.016 

0.039 

2 

0.001 

0.003 

0.005 

0.016 

0.039 

3 

0.002 

0.006 

0.01 

0.032 

0.078 

4 

0.003 

0.009 

0.015 

0.048 

0.117 

The  value  of  the  equivalent  network  capacitance  is. 


cx,r  •  «  vx,r  f,rais 


The  values  of  the  "conductances"  representing  conduction  between  elements 
are  presented  in  Table  III  in  terms  of  the  element  numbers  as  shown  in  Fig 
ure  5. 


For  the  circuit  response  calculation  the  linear  series  conduc¬ 
tances  can  be  combined  to  reduce  the  number  of  elements.  The  equivalent 
circuit  used  is  shown  in  Figure  6  with  the  parameter  values  used  listed  in 
Table  IV. 


Two  calculations  were  made  with  the  use  of  this  model  to  get  some 
insight  into  the  contribution  of  two-dimensional  carrier  flow:  one  case  was 
with  the  full  network  and  the  other  was  a  one-dimensional  network  with  the 
lateral  diffusion  conductances  were  removed.  In  both  cases  the  carrier 
generation  terms  were  assumed  to  be  a  11  ns  pulse.  This,  of  course.  Is  too 
long  to  be  a  realistic  single  particle  ionization  but  the  relative  one-  and 
two-dimensional  effects  should  still  be  illustrated. 


The  calculated  transient  Junction  current  is  shown  In  Figure  7  for 
the  one-  and  two-dimensional  cases.  It  Is  seen  that  the  variation  between 
the  two  Is  small  (-10%)  and  the  one-dimensional  approximation  used  In  many 
analyses  seems  reasonable. 
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Table  III.  Luaped  Model  Conductances  In  Units  of  qirL  (nhos). 


Gsi 

3  0.2 

G52 

=  0.2 

G53 

3  0.0667 

G54 

3  0.0667 

Gss 

3  0.04 

G56 

3  0.04 

Gs/ 

3  0.025 

Gss 

=  0.025 

G59 

3  0.0154 

G60 

3  0.0154 

G110 

3  0.2 

G120 

3  0.2 

Gl  30 

3  0.4 

Gl40 

3  0.6 

G6i  3  0.6 
G62  3  0.6 
Gg3  *  0.2 
G$4  3  0.2 
Ges  -  0.12 
Ggs  3  0.12 
Gs/  3  0.075 
Gee  3  0.075 
Ggg  *  0.046 
G/o  3  0.046 


Gioi  *  0.4 
Gin  *  0.4 
G121  3  0.4 
G131  *  0.8 

Gim  3  1.2 


G/l  3  1.0 
G/2  3  1.0 
G/3  3  0.33 
G/4  3  0.33 
G/s  3  0.20 
G/6  3  0.20 
G//  3  0.125 
G/s  3  0.125 
G/9  =  0.077 
Gso  3  0.077 


G102  =  0.4 
Gi 12  3  0.4 
Gl22  3  0.4 
6l32  3  0.8 
Gl42  3  1.2 


Gsi  3  3.2 
G82  3  3.2 
Gs3  3  1.067 
G84  3  1.067 
G85  3  0.64 
G86  3  0.64 
Gs/  3  0.40 
Ggs  3  0.40 
G89  3  0.246 
G90  3  0.246 


G103  3  0.4 
Gi 1 3  3  0.4 
Gl23  3  0.4 

Gi 3  3  3  0.8 

G143  3  1.2 


G91  3  7.8 
G92  3  7.8 
G93  3  2.6 
G94  =  2.6 
G95  3  1-56 
Ggs  3  1.56 
Gg/  3  0.975 
Ggs  3  0.975 
G99  3  0.6 
G100  3  0.6 


**&€*#**  ♦  „V'»w*S6Vl 


Table  IV.  Circuit  Description  Single  Event  2-d  Photoresponse. 


GS1  =  0.04  mhos 

Gsi i  s  0.01  mhos 
G512  *  0.005  mhos 
G513  s  0.0031  mhos 
G514  *  0.002  mhos 

G100  3  0.04  mhos 
Olio  *  0*04  mhos 
G120  =  0.04  mhos 
G130  *  0.08  mhos 
Gmo  *  0.12  mhos 

G61  *  0.12  mhos 
G611  *  0.03  mhos 
G6i2  a  0.015  mhos 
Ggi3  =  0.0092  mhos 
G611f  *  0.0058  mhos 

Gl01  *  0.08  mhos 

Gm  =  0.08  mhos 
G121  *  0.08  mhos 

G131  =  0.16  mhos 
G141  »  0.24  mhos 

Gn  »  0.2  mhos 
O/11  *  0.05  mhos 
O/12  *  0.025  mhos 
O/13  *  0.0154  mhos 


C00  3  2.0  x  10-10  F  =  C01  =  C02 

loo  =  1  X  10-3  A  .  f(t)  -  Iqi  3  I02 

Gqq  3  0.0002  mhos  3  Gq^  =  0g2 

C03  =  4  X  10-10  F 

I03  =  2  X  10-3  A  •  f(t) 

G0 3  =  0.0004  mhos 
C0„  -  6.0  «  10-10  F 
I„  =  3  «  10-3  A  .  fit) 

G04  =  0.0006  mhos 


C10  3  6  x  10-1°  F  3  Cn  3  C12 

G10  =  0.0006  mhos  3  G11  3  G12 

Cl3  "  1.2  *  109  F 

S13  =  0.0012  mhos 

C14  3  1.8  X  10-9  F 

Gm  3  0.0018  mhos 

C2o  3  1  *  10-9  F  3  C21  3  C22 

G20  3  0.001  mhos  3  G21  3  G22 

C23  3  2  x  10-9  F 

G23  3  0.002  mhos 

C2tt  3  3  x  10-9  F 

G24  3  0.003  mhos 
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Table  IV.  Continued. 


G81  =0.64  mhos 

C30  = 

3.2  X  10-9  F  =  C31  =  C32 

G8ii  =  0.16  mhos 

S30  = 

0.0032  mhos  =  G31  =  G32 

G8i2  =  0*08  mhos 

C33  = 

6.4  x  10-9  F 

Gsi 3  =  0.049  mhos 

G33  = 

0.0064  mhos 

G8m  =  0.03  mhos 

O34 

G34  = 

9.6  x  10-9  F 

0.0096  mhos 

G102  *  0.08  mhos 

C40  a 

7.8  x  10-9  F  =  C4I  =  C42 

G112  =0.08  mhos 

G40  = 

0.078  mhos  =  G4l  =  G42 

G122  a  0.08  mhos 

C43  = 

1.56  x  10-8  F 

G132  =  0.16  mhos 

G4  3  = 

0.0156  mhos 

G142  =  0.24  mhos 

C44  = 

2.34  x  10-8  F 

Gg  i  =  1.56  mhos 

&103 

=0.08  mhos 

G9U  =  0.39  mhos 

G1 13 

=  0.08  mhos 

G912  *  0.195  mhos 

0l23 

=  0.08  mhos 

G913  =  0.12  mhos 

G1  33 

=0.16  mhos 

Ggm  =  0.074  mhos 

0l43 

=  0.24  mhos 

Figure  7.  Calculated  transient  junction  current 
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